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The production of J/t/j mesons in photon photon collisions allows for a test of factorization 
in color-octet processes as predicted by NRQCD and observed at the Tevatron. We calculate 
the cross sections for J/tp production with direct and resolved photons, including also the 
feed-down from Xd and decays. Our NRQCD predictions are nicely confirmed by recent 
data from the DELPHI collaboration at CERN LEP2. 



1 Motivation 

One of the most interesting features of quantum chromo dynamics (QCD) is the apparent discrep- 
ancy between the fundamental role of color as a conserved quantum number carried by quarks 
and gluons and the non-observation of color in physical hadronic states. As a non-perturbative 
phenomenon, the confinement mechanism responsible for the color-singlet nature of hadrons is 
still not fully understood. One can hope to gain more insight by studying the relatively simple 
bound states of heavy charm and bottom quarks, such as J/ip or T mesons. 

Ever since QCD emerged as the fundamental theory of the strong interaction, models have 
been designed to explain t he com pensation of color in heavy quarkoni a. W ell-known examples 
are the Color Evaporation li^^l^ and Hard Comover Scattering Models^^^^, which assume that 
color is restored by soft and hard gluon exchange with the underlying event, respectively, and the 
Color Singlet Model (CSM)EIII, which requires the heavy quark pair to be produced in a color 
singlet state. The CSM leads to strong selection rules and connects the production cross section 
with the quarkonium wave function. Today, a rigorous effective field theory exists in the form 
of non- relativistic QCD (NRQCD) 13, which utilizes a double expansion in the strong coupling 
constant and the relative quark-antiquark velocity v in order to factorize the hard production 
from the soft binding process^. Left-over singularities in the CSM can be removed systematically 
into non-perturbative color-octet operator matrix elements (OMEs), and their numerical values 




Figure 1: Feynman diagrams pertinent to the partonic subprocess 77 J/V'7- 

can be fitted to describe the hadroproduction cross sections observed at the Tevatron. However, 
it is necessary to demonstrate the universahty of these OMEs in other production processes, 
such as 77 collisions. Since the theoretical uncertainties from scale variations in leading order 
(LO) of Us are quite substantial, it is furthermore important to include also contributions from 
next-to-leading order (NLO) virtual loop and real emission processes. In this paper, we report 
on recent progress along these lines. 

2 Real corrections to J production in direct 77 collisions 

In direct 77 collisions, J/ij: mesons with mass M = 2mc and finite transverse momentum px are 
produced in association with either a photon (see Fig. ^ or a jet, where the final state photon 
in Fig. ^has to be replaced by a gluon. In the first case, the physical color-singlet state {^S^^^ in 
spectroscopic notation) is produced directly, whereas in the second case the intermediate color- 
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octet state ( 5| ) transforms into the physical J/ip through soft gluon emission. Although the 
strong coupling of the hard gluon enhances the second process, it is strongly suppressed by the 
relevant color-octet OME, which is subleading in v. In NLO of as, an unresolved dijet system 
with invariant mass sjj, originating from two gluons or a light quark-antiquark pair, allows for 

the production of intermediate ^Sq'\ '^Si'\ and ^Pj"^ states. At small pr, virtual loop corrections 
have to be included to cancel the unphysical dependence on the cut-off Sjj > M^, but at large 
Pt these contributions become unimportant. As can be observed clearly in Fig. [2 the cut-off 

ee ^ J/V|/ + X at LEP2 (Bremsstr.) 




Figure 2: Transverse-momentum distribution da/dpr of 77 3 j^^X via bremsstrahlung at CERN LEP2. The 
sum of the LO contributions for X = ^ and X = 7 is compared with the 2 — > 3 part of the NLO contribution for 
dijet mass Sjj > and MV20. For comparison, also the '^5'[*' -channel contributions to the latter are shown. 

dependence is reduced and the corrections from real particle emission become larg At small 
Pt, one does, however, not expect large radiative corrections^ as has also been shown in a full 



NLO calculation for direct color-singlet J/V' photoproductioni 



3 Comparison of J/ifj production in complete 77 collisions with CERN LEP2 data 



J/^jj production in 77 collisions proceeds not only through direct interactions of photons with 
charm quarks, but also through resolved processes, where one or two of the photons emit quarks 
and gluons which then interact with the heavy quarks. In the small-pT range recently studied by 
the DELPHI experiment at CERN LEP2lI^, J/ ip production is dominated by the single-resolved 
process 'jg — > J/ipg, where one initial state photon and the final state photon in Fig. ^ have 
to be replaced with gluons. In addition, processes with non-abelian coupling of the gluons and 
processes involving light quarks have to be taken into account using the relevant color-octet 
OMEs. Furthermore, contributions from Xcj and ip' decays have to be included, whereas those 
from B meson decays are suppressed by the small branching fraction. 

In our numerical analysis, we use nic = (1.5 ± 0.1) GeV, a = 1/137.036, and the LO 
formula for a^r^\fJ-) with Uf = 3 active quark flavors. As for the photon PDFs, we use the 
LO set from Gliick, Reya, and Schienbein (CRS)*^, which is the only available one that is 
implemented in the fixed-flavor-number scheme, with Uf = 3. We choose the renormalization 



and factorization scales to be = i^mx and /x/ = ifrnx, respectively, where mx = \J M'^ + Pr 
is the transverse mass of the J/ip meson, and independently vary the scale parameters and 
between 1/2 and 2 about the default value 1. As for the J/V'i Xcj-, and ip' OMEs, we adopt the 
set determined in Ref.^1 by fitting the Teva tron data using the LO proton PDFs from Martin, 
Roberts, Stirling, and Thorne (MRST98LO)l^ as our default and the one referring to the LO 
proton PDFs from the CTEQ Collaboration (CTEQ5L)II3 for comparison (see Table I in Ref. 
^1). In the first (second) case, we employ Aq^j) = 204 MeV (224 MeV), which corresponds to 
AgJ,j3 = 174 MeV^ (192 Me\^, so as to conform with the fit . Incidentally, the GRS 
photon PDFs are also implemented with Aqqj-, = 204 MeV^^. In the cases ip = J/ip^ip', the 
fit results for [O^'^sf^]) and {0^fpf\) are strongly correlated, and one is only sensitive to the 
linear combination Mf = {0't'[^sf]) + with an appropriate value of r. Since the 

cross section is sensitive to a hnear combination of (©^[^^f^]) and {0^[^pf^]) that differs from 
Mf, we write {0'*[^sf^]) = nMf and (O^ppj^^]) = (1 - k) {ml/r)Mf and vary k between 
and 1 about the default value 1/2. The J-dependent OMEs (O^ppj^']), {O^-J [^P^j\) , and 
satisfy multiplicity relations, which follow to leading order in v from heavy-quark 
spin symmetry. In order to estimate the theoretical uncertainties in our predictions, we vary 
the unphysical parameters .^j, and k as indicated above, take into account the experimental 
errors on m-c, the decay branching fractions, and the default OMEs, and switch from our default 
OME set to the CTEQ5L one, properly adjusting Aq^j-,. We then combine the individual shifts 
in quadrature, allowing for the upper and lower half-errors to be different. 

In Fig.|31we confront the p^^ distribution of e^e~ — > e^e~ J/ip-\-X measured by DELPhJI^I 
with our NRQCD and CSM predictions. The solid lines and shaded bands represent the central 
results, evaluated with our default settings, and their uncertainties. The DELPHI data clearly 
favors the NRQCD prediction, while it significantly overshoots the CSM one^^. 



4 Conclusion 

In summary, the production of J/tp mesons in 77 collisions allows for a test of factorization in 
color-octet processes as predicted by NRQCD and observed at the Tevatron. We have calculated 
the contributions from direct and resolved photons to J/ ip mesons produced directly or from Xcj 
and ip' decays. Our NRQCD predictions are nicely confirmed by recent data from the DELPHI 
collaboration at CERN LEP2, whereas the CSM prediction is clearly disfavored. 
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Figure 3: The cross section da/dp^ of e"'"e ^ e^e J/ip + X measured by DELPHl'— 'as a function of is 
compared with the theoretical predictions of NRQCD and the CSM. The soli d and dashed lines rep rese nt the 
central predictions obtained with the ME sets referring to the MRST98LoE^ (default) and CTEQSlEU PDFs, 
respectively, while the shaded bands indicate the theoretical uncertainties on the default predictions. The arrows 
indicate the NRQCD prediction for pt = and its ^ , ^S"!,*' , ^5*^ , and ^Pj*' components. 



Acknowledgment 

The author thanks the organizers of the QCD Moriond 2003 conference for the kind invitation 
and B.A. Kniehl, L. Mihaila, and M. Steinhauser for a fruitful cohaboration. 

References 

1. H. Fritzsch, Phys. Lett. B 67 (1977) 217. 

2. F. Halzen, Phys. Lett. B 69 (1977) 105. 

3. M. Gliick and E. Reya, Phys. Lett. B 79 (1978) 453. 

4. P. Hoyer and S. Peigne, Phys. Rev. D 59 (1999) 034011. 

5. N. Marchal, S. Peigne and P. Hoyer, Phys. Rev. D 62 (2000) 114001. 

6. E. L. Berger and D. L. Jones, Phys. Rev. D 23 (1981) 1521. 

7. R. Baier and R. Riickl, Phys. Lett. B 102 (1981) 364. 

8. W. E. Caswell and G. P. Lepage, Phys. Lett. B 167 (1986) 437. 

9. G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. D 51 (1995) 1125 [Erratum-ibid. 
D 55 (1997) 5853]. 

10. M. Klasen, B. A. Kniehl, L. Mihaila and M. Steinhauser, Nucl. Phys. B 609 (2001) 518. 

11. M. Kramer, J. Zunft, J. Steegborn and P. M. Zerwas, Phys. Lett. B 348 (1995) 657. 



12. S. To dor ova-Nova, hep-ph/0112050; M. Chapkin, private communication. 

13. S. Frixione, M. L. Mangano, P. Nason and G. Ridolfi, Phys. Lett. B 319 (1993) 339. 

14. M. Gliick, E. Reya and I. Schienbein, Phys. Rev. D 60 (1999) 054019 [Erratum-ibid. D 62 
(2000) 019902]. 

15. E. Braaten, B. Kniehl and J. Lee, Phys. Rev. D 62 (2000) 094005. 

16. A. Martin, R. Roberts, W. Stirling and R. Thorne, Eur. Phys. J. C 4 (1998) 463. 

17. H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. G 12 (2000) 375. 

18. M. Klasen, B. Kniehl, L. Mihaila and M. Steinhauser, Phys. Rev. Lett. 89 (2002) 032001. 



